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INTRODUCTION 

A number  of  problems  are  associated  with  describing  the  propagation 
of  ultrasonic  waves  in  laminated  materials.  In  general  such  materials  are 
both  inhomogeneous  and  anisotropic;  therefore,  there  exist  a number  of  modes 
of  wave  propagation  associated  with  different  directions  of  propagation. 
Attempts  to  treat  these  materials  as  monolithic  materials  with  effective^ 
elastic  moduli  provide  reasonable  descriptions  of  quasi-static  behavior;' 
however,  when  one  attempts  to  describe  transient  phenomena  with  such  a 
model,  it  is  found  that  no  choice  of  effective  parameters  is  adequate  to 
describe  the  high  frequency  behavior  of  propagating  waves  as  influenced  by 
tht:  periodic  nature  of  the  composite.  Recently  a number  of  investigators^’  ^ 
have  approached  the  problem  of  waves  propagating  In  a composite  medium 
utilizing  Fluoquet's  Theorem  for  describing  wave  propagation  in 
infinitely  periodic  structures.  Most  of  this  work  on  infinite  laminates 
has  been  carried  out  in  order  to  predict  wave  speeds  associated  with  the 
effect  of  impact  on  laminates.  These  results  have  been  promising;  however, 
it  is  not  clear  to  what  extent  they  can  be  applied  to  finite  periodic 
structures. 


The  problem  of  wave  propagation  in  finite  laminates  has  been 
considered  by  a number  of  Investigators  in  the  field  of  nondestructive 
testing,^*  5 who  have  been  investigating  the  effect  of  laminates  on  the 
propagation  of  ultrasonic  pulses,  in  order  to  develop  nondestructive  testing 
techniques  for  probing  the  structure  of  these  laminates. 


A number  of  investigators  have  observed  periodic  fluctuations  in 
the  frequency  transmission  and  reflection  spectra  associated  with  laminated 
structures.  This  periodicity  In  the  frequency  spectrum  has  been  explained” 
as  resulting  from  interference  of  coherent  reflections  of  ultrasonic  pulses 
from  two  separated  Interfaces  within  a sample. 


It  will  be  shown  that  this  two  interface  model,  while  phenomeno- 
logically correct  In  some  cases,  fails  in  general  to  accurately  describe 
wave  propagation  in  laminates  in  which  several  interfaces  exist.  A 
modified  model  will  be  derived  here  v/hich  is  capable  of  dealing  with  the 
general  case  of  finite  laminates  and  their  ultrasonic  frequency  spectra. 

Also,  the  general  interpretation  of  frequency  spectra  in  laminates  will  be 
discussed  with  emphasis  on  applications  of  spectrum  analysis  to  nondestructive 
testing  of  laminates  and  composites. 
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THEORY 


A.  Acoustic  Waves  Incident  on  a Boundary 


It  is  a well  known  result^  that  In  the  absence  of  dispersion  the 
transmission  and  reflection  amplitude  coefficients  T and  R of  acoustic 
plane  waves  incident  normally  from  a medium  with  acoustic  impedance  Z] 
onto  a medium  of  acoustic  irr^edance  Z2  are  given  by: 


(la)  R 


12 


-Z 


(lb)  T 


12 


2Z. 


Z2  -.z^ 


where  the  frequency  independent  Z for  each  medium  is  given  by  the  product 
of  the  density  p and  the  speed  of  sound  C. 

When  more  than  two  media  are  involved,  interference  effects  produce 
frequency  dependent  transmission  and  reflection  coefficients  for  waves 
traveling  through  media  of  finite  thickness. 

In  general,  for  a given  frequency,  steady  state  solutions  exist  which 
can  satisfy  the  boundary  conditions  for  wave  propagation  through  arbitrary 
laminar  media;  these  solutions  may  be  obtained  through  the  solution  of  2n-l 
Pnear  simultaneous  equations.  The  result  of  solving  these  equations  is 
to  generate  steady  state  frequency  dependent  transmission  and  reflection 
coefficients.  T((o)  and  R (u) . 


Because  of  the  linear  nature  of  the  governing  equations,  implicit 
solutions  to  transient  problems  can  be  constructed  from  these  transmission 
and  reflection  coefficients  using  a Fourier  integral  technique.  For 
example  If  a number  of  layered  media  are  impacted  by  a finite  v;ave  train 
of  amplitude  A(t)  and  Fourier  transform  A(u) , the  transmitted  frequency 
and  time  depervdent  wave  forms  will  be  given  by: 

(t)  “ ^e  * T(u)  A(w)  duj 


Ay  (w)  " T (u)  A(<i>) 

and  the  spectrum  A^^  ^2  (w)  of  any  time  slice  (t2  - tj)  of  that  transmitted 
wave  train  Is  given  by: 


-1 

2 


\lt2  ^ *** 

J t2 


tj  >tj  > t=o 


where  t ■ o is  the  time  of  arrival  of  the  wave  train  at  the  first  interface. 
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Both  A (»  ) (the  spectrum  of  the  total  transmitted  wave  form) 
and  A (u,)  ^ (the  spectrum  of  the  gated  wave  form),  can  be 

experimentally  determined  from  ultrasonic  measurements. 

In  this  paper  we  will  develop  expressions  for  the  quantities 
T {<t  ) and  R (to  ) by  directly  applying  the  principle  of  superposition 
to  determine  the  response  of  a layered  medium  to  a transient  wave  train. 

This  approach  will  generate  a summabte  infinite  series  for 
R ( <0  ) and  T ( u );  in  addition  to  expressions  for  A^  from 

finite  partial  sums.  An  iteration  scheme  will  be  i 2 developed 
for  calculating  exact  values  of  R (u  ) and  T (<«  ) for  layered  media. 
These  values  will  be  built  up  from  the  results  for  a single  layer,  and 
each  successive  iteration  will  add  an  additional  layer  to  the  medium. 

Thus  in  computing  T ( o> ) and  R ( u)  for  n layered  media,  results  for 
n-1,  n-2,  n-3,  etc.  layers  will  be  simultaneously  generated. 

B.  Wave  Propagation  in  a Single  Layer 

Consider  the  getxnetry  of  Figure  1.  Let  waves  Incident  from  the 
left  upon  medium  2 have  transmission  and  reflection  coefficients  Tj2  and 
Rj2  and  those  incident  from  the  right  T£|  and  R2^.  The  coefficients  for 
tne  boundary  between  medium  S and  3 are  analogously  defined* 

A portion  of  the  pulse  A(t)  lr'-*'*^nt  from  the  left  onto  medium  2 
will  be  transmitted  with  att^litude  T|.  ' and,  i’’>on  emerging  in  medium  3» 

will  produce  a disturbance  of  the  form 


'’’23  '’’12  A(t  - £ ) 

in  addition  to  the  portion  of  the  wave  transmitted  into  medium  2, 
there  will  be  a nund>er  of  pulses  onerging  which  have  undergone  multiple 
reflections  within  this  medium.  The  total  of  all  of  these  waves  is 
given  by: 

^3  d A (t-  1 ) + T,2  R23  Rj,  T23  A (t  - M ) 

(2)  a’  (t)  • T,2  T23  Y.  '^21  )*^''<«‘<2N+0d) 

where  c • 

Similarly  the  total  wave  reflected  back  into  medium  1 is  gi^n  by: 

' d “1  ^ 

(3)  R I (t)  2 A{t)  + T,2  T^j  R2,  £ (R23  ^21^  ^ 

^ c 

N-1 
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Taking  the  Fourier  transforms  of  expressions  (2)  and  (3)  and 
utilizing  the  Identity 


Jf  (t  -At)  e dt  - e Jf(t) 


i«t  . 
e dt 


M have: 

W A'  { to)  - Ft  T e 
13  L 12  23 


-lud'/c^ 


(5)  R;<»)  - (R,J  ♦ e 


l/c«  u -2lu«<l/c  1 

T,  (%,  %,)  ‘ ■!  *(“) 

(4^0  23  21 

-2i«d/c  ^ , -2Iojd/c  N 

^ («23  "2,  e ) A(„)  ) 


I * 

where  A (u)  Is  the  Fourier  transfdrm  of  A^^  (t) , and  R ^ (to)  Is 

the  Fourier  transform  of  R^  (t)»  etc. 

The  sum  appearing  In  expressions  A and  5 can  be  readily  evaluated  for 
N either  finite  or  infinite  allowing  the  definition  of  generalized 
frequency  dependent  reflection  and  transmission  coefficients  r(u)  and  t{u) 
for  the  layered  medium.  For  R Infinite  we  have  the  expression 


(6)  t(u)  i a'  (u)  « T^2  ® 


-iud/c 


1 - R,,  h]  e 


(7)  r{o>) 


= (u) 

A (u) 


T,_  1 R e 
12  21  23_ 


' - "23  "21 


■2 i ud/c 


— + R,, 

-2io)d/c  *2 


C.  Attenuation  in  a Layer 

The  use  of  complex  notation  allows  attenuation  effects  to  be  introduced 
into  equations  (6)  and  (7).  This  is  done  most  easily  by  allowing  d to  be 
complex. 

d-  d (1  -lo) 

where  d is  a real  distance  and  ct  becomes  a linear  attenuation  coefficient, 
corresponding  to  a wave  whose  amplitude  is  exponentially  damped  with 
distance. 

0.  Generalized  Reflection  and  Transmission  Coefficients 


Clearly,  for  the  case  of  layered  media  of  finite  thickness  a single 
reflection  coefficient  does  not  exist  for  the  time  dependent  representation 
Sf  the  wave  A (t);  howerer,  we  can  define  operators  5 and  f which  when 
(grating  on  A<t)  g«ierate  the  expressions  for  tl«  tiw  dependent  reflected 
mad  transmitted  waves  respectively.  The  appreprlete  definition  for  rhese 
opereticms  le: 
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(8)  R A(t)  - 

1/2  r -1/2 

- (2H>  l(t(u)  exp  (-Ittt)  (2n) 

* “l-«o 

(9)  TA(t)-  ^ ^ (tM  ^(A  (t*))) 

1/2  f -1/2 

(2H)  I (t(w)  exp  (-Iwt)  <2n) 


a 

1 


A(t')  expClitt*)  dt')dt» 


a 

f 


A(t')  e3q>  (i6Jt)dt')dB 


wkmf  ^ and  ^ are  the  Fourier  transform  and  its  inverse.  1!he  use 
of  this  operator  notation  permits  a layer  to  be  represented  as  a boundary 
or  virtual  interface  having  these  operators  as  their  reflection  and 
transsLasim  coefficients  (see  Figure  2). 


E.  Wave  Propagation  in  Multi-Layered  Media 

Consider  the  layered  media  In  Figure  2.  This  is  essentially  the 
same  as  Figure  1 except  that  the  boundary  between  media  2 and  3 is  now 
replaced  by  a boundary  medium  with  tranmisslon  and  reflection  operators 


®23 


and 


23 


Rewriting  expressions  (2)  and  (3)  in  the  form  of  a series  and 


renlaclns  T 


23 


and  by  operators  we  have* 


(10)  A* 


13 


(11)  R’ 
1 


(t) 

A(t)  + T 


23 


oo 

2. 


N 

(R  R ) A(t-(2iH-l)d/c  ) . 
23  21 


12 


T R 
12  21  23 


00  N 

V (R  R ) A(t-2Nd/c) 
23  21 


, It  is  readily  verified  by  recourse  to  their  definitions  that  rhe 
transmission  and  reflection  operators  coa»ite  and  that  the  Fourier  transform 
of  a number  of  operators  acting  on  a function  is  just  the  product  of  the 
corresponding  frequency  reflection  and  transmission  coefficients  times 
the  Fourier  crsasforas  of  the  function. 


T R Aft)  - R T A(t) 

N N N N 

and  ^(R  T A(t))-  r (u)  t («)  ^ (A(t)) 

Applying  these  results  to  take  the  transform  of  expression  (10)  «id  (11) 
ve  ha'». 


(12)  T'  (»)  - A’  (u)  / A Cm)  “ T tCa)  (r(6>)  R ) exp  (-lca(2!«-l)d/c) 
13  12  ^ 21 

" N-0  N 

(13  R’  («)  ' R (b)  + T T r (a)  > (f  (a)  R ) exp  (-2i!feal/c) 

12  10  91  91 


12  21 
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Again  summing  the  series  we  have, 

, . ' -iwd/c 

(12a)  T (a.)  - Ti2  t (a,)  e 

1 “ r (ta)  R2t  e IwSTc 

(Ua)  r'  (o,)  - R,2  + T^2  r(o)) 

^ “ *‘l  M *^21 


-2io)d/c 

e 


e “2iu)d/c 


Successive  application  of  these  formulae  yields  a simple  iterative  proce- 
dure for  tabulating  transmission  and  reflection  coefficients  for  media  having  an 
arbitrary  number  of  layers  of  varying  thickness  and  acoustic  parameters. 

Normal  incidence  reflection  and  transmission  coefficients  for 
several  finite  arrays  have  been  calculated  using  formulae  (la),  (lb),  (12a) 
and  (12b)  and  febe  results  of  these  calculations  are  discussed  below. 


F.  Calculations  of  Reflection  and  Transmission  Spectra 


1.  Single  Laminae  in  an  Ambient  Medium 


The  values  of  the  functions  log  joL^(w)l  loKCw)!  ® 

0.30"  steel  plate  In  water  are  plotted  in  Figures  3a  and  3b  respectively. 
The  values  of  these  coefficients  were  computed  directly  from  expressions 
(12a)  and  (13a)  using  a Fortran  computer  program  "Layer"  written  for  this 
purpose.  A computer  controlled  plotter  was  used  to  reproduce  the  plots 
directly  from  computed  values  with  a resolution  of  100  points  per  inch  on 
the  frequency  axis.  This  resolution,  although  acceptable  for  plotting 


10 
log 


t((o) 

r(u) 


is  inadequate  to  reveal  the  sharp  resonances  In  the  function 
which  is  in  fact  periodic.  The  existence  of  these  sharp 


1A  I * \w/  I I VI*  99  III  low  poiiv/uiv*  I lie  e/vi9W^Mve  9norp 

periMically  spaced  resonances  termed  thickness  resonances  (period  f*  c/2d  ) 
can  be  exploited  for  highly  accurate  sound  velocity  measurements. 


2.  Multi-layered  Periodic  Media 

In  addition  to  reproducing  the  known  analytical  results  for  trans- 
mission and  reflection  from  single  layers,  "LAYER"  has  been  applied  to 
calculating  1o^q|  t(o)),|  and  log  jQ|r(u)!  for  finite  periodic  layered  arrays. 
The  results  of  «ych  calculations  are  seen  In  Figures  3c  through  3g.  which 
plot  log  IQ  |t(w)|for  stacks  consisting  of  ! to  5 glass  plates  1.22mm  thick 
spaced  .Imm  apart  in  an  ambient  medium  of  water.  It  is  interesting  to  note 
that  when  more  than  a single  plate  is  present  log  ^q|  |t(u)|  is  no  longer  a 
simply  periodic  function. 


The  shallow  tightly  spaced  minima  seen  in  Figures  3d  through  3g 
correspond  roughly  to  the  thickness  reasonance  of  the  entire  array  and  they 
increase  In  nuBiber  as  plates  are  added  to  the  array.  However,  the  broader 
resonances  remain  constimt  in  number  but  rapidly  deepen  as  layers  are  added, 
producing  frequency  bands  having  increasingly  high  acoustic  attenuation. 
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3.  Infinite  Periodic  Media 

It  can  be  shown  that  In  the  limit  of  a periodic  medium  with  an 
Infinite  number  of  layers,  there  will  be  frequencies  for  which  no  traveling 
plane  wave  solutions  exist  inside  the  medium.  This  is  a limiting  case  of 
the  behavior  seen  In  the  transmission  curves  in  Figures  3c  through  3g. 
Robinson  and  Leppelmeler  (Journal  Applied  Mechanics,  March  197^«  p*  89) 
have  studied  the  case  of  wave  propagation  In  a two  component  layered 
periodic  medium  and  have  given  an  expression  for  the  wave  number  k of 
longitudinal  traveling  waves  propagating  in  such  a medium: 


(Ik)  cos  kl  ■ coe  («  d^)  cos  (u»  d2)  * 1,  (pj  + p C^)  slntw  d (“  ^2^ 

^2  ^2^2  ® 


'1 


Cl  and 


are  the  layer  thicknesses  for  media  having  sound  velocities 


where  d.  and  d2 

C2  respectively  and  !<■  d^  + do 


Uj  mu  i#2  reapc&Liwiy  aiiu  3«<*  u]  -r  Whenever  p.  y p2  C2  in  expression 
l4,  there  exist  values  of  k for  which  {cos  Kll>  1,'  This  implies  a complex 
value  of  k and  frequency  regions  corresponding  to  such  wave  numbers  are 
referred  to  as  forbidden  bands,  because  there  exist  no  traveling  waves 
propagating  at  these  frequencies.  Using  expression  (Ik),  the  first  four 
forbidden  frequency  bands  have  been  calculated  for  the  layered  glass  array 
of  Figures  3c  - 3g  and  are  indicated  by  arrows  In  Figure  3g*  It  is  seen 
that  these  forbidden  bands  correspond  to  broad  transmission  minima  which  are 
seen  to  deepen  rapidly  as  the  number  of  layers  in  the  finite  array  is 
increased. 


Experiment 

Acoustic  reflections  and  transmission  coefficients  were  measured  for 
a number  of  materials  using  the  apparatus  pictured  schematically  in  Figure  k. 
Broadband  longitudinal  ultrasonic  stress  waves  a r»  projected  onto  a sample  by 
a pulsed  piezoelectric  trandsucer  (A).  The  projected  pulse  Interacts  with 
the  sample  producing  transmitted  and  reflected  wave  trains.  The  transmitted 
wave  train  Is  detected  by  a second  transducer  (B)  placed  behind  the  specimen, 
while  the  reflected  wave  Is  detected  by  the  projecting  transducer. 

Either  the  transmitted  or  the  reflected  Dulse  mav  be  selected  for 
amplification  and  gated  electronically  to  obtain  the  desired  portion  or  time 
slice  of  the  signal.  This  gated  signal  is  processed  by  a spectrum  analyzer, 
which  generiUi  a plot  of  the  amplitude  cf  the  gated  signal  as  a function  of 
frequency,  (l.e.  I\it2 
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RESULTS  AND  DISCUSSION 

A.  Transmission  and  Reflection  Through  a Single  Finite  Medium 

Figure  5 shows  the  results  of  an  experimental  measurement  of 
the  frequency  dependent  transmission  coefficient,  t(w),  for  a 0.3(f*  steel 
plate.  The  projecting  transducer  In  this  measurement  was  a Panametrics 
VIP  15  MHz  broadband  device  with  a 1/4“  diameter  piezoelectric  element;  the 
receiver  was  a 25  MHz  Automation  Industries  transducer  with  a 1/2“  diameter 
element.  The  distance  betwwen  the  sample  and  the  receiving  transducer  was 
2“.  The  gate  was  set  to  receive  all  of  the  transmitted  pulse  for  which  the 
amplitude  was  greater  than  the  ambient  noise. 

The  upper  reference  curve  In  Figure  5 Is  a plot  of  the 
logarithm  base  10  of  the  amplitude  of  the  signal  oroduced  at  the  receiving 
transducer  vs.  frequency  (with  no  sample  presend«  The  lower  curve  Is  a plot 
of  the  logarithm  of  the  transmitted  amplitude  as  a function  of  frequency. 

Tlie  spectrum  has  a bandwith  of  1 KHz  filtered  to  remove  most  of  the  rapid 
oscillations  associated  with  the  pulse  repetition  rate  of  the  pulser 
receiver.  Using  this  type  of  plot,  the  logarithm  of  the  transmission 
coefficient  can  be  measured  as  the  difference  In  the  ordinates  of  the  two 
curves. 

B.  Velocity  of  Sound  Measurement 

The  theoretical  values  #Or  the  transmission  coefficient  of  a 
0.3“  steel  plate  computed  from  the  absolute  value  of  expression  (6)  are 
plotted  In  Figure  3a.  The  function  Is  ,jer!odIc  with  period 

(15a)  Af„-C 
N 2d 

Frequencies  for  which  (15a)  is  satisfied  are  often  referred  to  as  thickness 
resonances.  Knowing  the  thickness  of  the  plate  this  formula  can  be  used  to 
measure  the  phase  velocity  of  sound  as  a function  of  frequency  from  experi- 
mental data  using  the  relationship. 

(I5b)  ' Ifl 

N 

For  steel  this  is  most  accurately  measured  using  reflectance  minima  for 
which  the  resonances  are  much  sharper  than  for  transmission  peaks  (see 
Figure  5)*  The  precision  of  this  technique  is  easily  better  than  one  part 
in  1000  for  steel;  hcv.*eyer,  the  yarlatlons  puhn«:hed  values  for  C make  it 
impossible  to  eonf I rm  these  values  to  better  than  2 significant  Figures. 
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C.  Multilayer  Media 

Figure  6 Is  a plot  of  the  transmission  spectrum  from  a layered 
array  of  three  glass  plates  .00I22M  thick  separated  by  .001  meters  of  water. 
The  entire  array  was  Immersed  In  water  during  the  measurement  and  the  system 
represents  a somewhat  Idealized  laminate.  Figure  3e  Is  a plot  of  the 
transmission  spectrum  for  this  laminate  as  computed  from  formulas  I2a  and 
1 3a.  The  agreement  between  the  theoretical  and  experimental  spectra 
indicates  that  the  linear  acoustic  approximation  used  Is  a reasonable  model 
for  wave  propagation  In  laminated  structures.  Interesting  features  of  the 
complicated  laminate  spectrum  are: 

1.  Transmission  peaks  do  not  recur  with  simple  periodicity 
as  they  do  with  monolithic  materials. 

2.  Although  peaks  approximately  corresponding  to  the  first 
two  thickness  resonance  modes  for  the  entire  array  do  appear,  the  third 
harmonic  Is  absent  and  the  next  maximum  appears  after  a large  frequency  gap. 

These  results  are  noteworthy  because  It  has  often  been  Indicated 
In  the  nondestructive  testing  literature  that  reflections  from  a lamlna-llke 
structure  within  a specimen  will  produce  a periodic  spectrum.  This  Is 
probably  true  only  when  othef  reflections  In  the  structure  are  small,  are 
of  an  essentially  incoherent  nature  or  can  be  Isolated  using  a gating 
technique. 

D.  Agreement  between  Theory  and  Experiment 

As  mentioned  above,  velocities  of  sound  measured  by  the  spectrum 
analysis  technique  have  been  found  to  be  in  good  agreement  with  published 
values.  In  addition,  when  published  values  for  pand  c are  used  In  equations 
(1),  (6)  and  (7)  to  compute  spectra  of  laminates,  we  are  able  to  compute 
positions  and  relative  magnitudes  of  extremely  detailed  features  of  these 
spectra.  (See  Figures  6 and  3e). 

A somewhat  more  conclusive  test  of  the  agreement  between  theory 
and  experiment  Is  the  prediction  and  -byolyte  values  of  trans- 

mission and  reflection  coefficients,  for  example,  tha  minimum  transmission 
coefficient  of  the  steel  layer  In  Figure  5 is  about  4dB  above  Its  predicted 
value  of  -24dB  and  the  maximum  transmission  coefficient,  which  should  reach 
OdB  varies  t 2dB  from  this  value. 

The  observed  discrepancy  between  theory  and  experiment  is 
probably  due  to  near  field  effects  of  the  transducer.  This  hypothesis  has 
been  partially  confirmed  by  experl«wents  which  have  shown  decreasing  agree- 
ment between  predicted  and  observed  frequency  spectra  as  frequency  is 
Increased.  Since  the  near  field  of  the  transducer  has  a radius  which  Is 
proportional  to  the  transducer  diameter  divided  by  the  square  of  the  wave- 
length, the  near  field  always  extends  further  at  hiah  freauencies. 

Efforts  to  work  entirely  In  the  far  field  of  the  transducer  have  been 
hampered  by  poor  signal  to  noise  ratios  and  by  the  physical  constraints  of 
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the  {mtnersion  tank.  Improved  agreement  between  theory  and  experiment  may 
require  the  employment  of  analytical  techniques  for  modeling  near  field 
transducer  effects. 

Additional  experimental  error  results  whenever  the  bandwidth 
of  the  spectrum  analyzer  Is  on  the  order  of  the  spacing  between  spectral 
features.  This  results  In  the  smearing  out  of  sharp  spectral  features  such 
as  the  transmission  peaks  in  Figure  5,  which  should  be  much  sharper  and 
should  reach  a higher  maximum  value. 

E.  Frequency  Spectra  from  Structural  Composite  Materials 

Figure  7 shows  a plot  of  the  frequency  spectra  from  three 
different  graph ite/epoxy  laminates.  For  each  of  these  spectra  reflectance 
minima  exist  which  correspond  roughly  to  thickness  resonances,  ({.e.  f=Mc/2d) . 
This  suggests  that  to  a first  approximation,  a graphite  epoxy  laminate  may 
be  modeled  as  a monolithic  material.  However,  it  is  seen  that  as  in  the 
case  of  the  layered  glass  laminate  the  portion  of  the  frequency  spectrum 
corresponding  to  the  thickness  resonance  of  a single  period  of  the  laminate 

behaves  anomalously  (I.e.  resonance  number  4 is  missinq  for  the  4 dIv  laminate, 
resonance  number  18  is  missing  for  the  18  ply  laminate).  Pursuing  the 
implications  of  modeling  the  composite  as  a monolithic  material,  formula 
(15b)  may  be  applied  to  compute  a dispersion  relation  c(f)  for  each 
frequency  at  which  a thickness  resonance  occurs.  Such  a dispei jlon  relation 
Is  plotted  in  Figure  8 for  an  18  ply  graph ite/epoxy  laminate.  Significant 
features  of  this  dispersion  relation  are, 

(1)  A general  increase  In  velocity  with  increasing  frequency, 

(a  phenomenon  expected  In  viscoelastic  materials). 


4 MHz. 


(2)  Finite  discontinuities  in  velocity  at  multiples  of 


F«  Forbidden  Frequency  Bands  in  Composites 

The  existence  of  discontinuities  in  the  dispersion  relations 
for  composite  materials  results  from  the  periodic  structure  of  the  composite 
and  may  be  predicted  from  equation  (14)  and  the  relation, 

c-f/2nk 

It  is  found  that  finite  discontinuities  !n  velocity  occur 
when  frequencies  of  high  attenuation  referred  to  as  fo.-*bldden  frequency  bands 
are  reached.  Alternately  may  say  that  these  discontinuities  occur  whenever 
the  wave  number  K is  a multiple  of  2n/i  where  i is  the  thickness  of  ons  period 
of  the  medium  . (See  Figure  7)>  Thus,  although  wave  propagation  does  not 
cease  in  the  forbidden  bands  for  finite  layered  laminates,  transmission  Is 

strongly  attenuated.  In  addition  if  experimental Iv  measured  dispersion  curves 
for  finite  media  are  examined,  it  is  found  that  regions  of  anomalous  disper'* 
sion  occur  In  the  vicinity  of  the  forbidden  bands. 
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It  is  interesting  to  note  that  the  graphite/epoxy  composites 
whose  spectra  are  plotted  In  Figure  7 are  not  1aid-up  in  a truely  periodic 
array  but  consists  of  0°,  90°  and  + 45°  laydrs.  Some  indication  of  the 
effects  of  ply  arrangment  can  be  seen  in  Figure  7 (Velocity  dispersion 
for  an  18  ply  composite).  In  addition  to  the  discontinuities  which  occur 
near  8 and  16  MHz  corresponding  to  resonances  of  a single  layer,  anomalies 
also  occur  near  4 and  12  MHz  which  correspond  to  resonances  of  lamina 
pairs.  Strictly  speaking  one  cannot  exaplain  such  resonances  In  terms  of 
the  simple  theory  presented  here,  since  for  normal  incidence  longitudinal 
waves  identical  angle  ply  laminae  of  different  orientation  should  behave 
identically.  It  appears  that  for  the  general  graphite/epoxy  laminate  a 
detailed  model  for  stress  wave  propagation  must  include  the  effects  of 
waves  being  scattered  !n  directions  other  than  parallel  to  the  incident 


G.  Applications  of  Spectrum  Analysis  in  Nondestructive  Testing 

The  reasonable  success  in  modelling  the  ultrasonic  frequency 
spectrum  of  layered  materials  and  the  empirically  confirmed  capability  of 
predicting  qualitatively  the  frequency  spectra  of  fiber  reinforced  composites 
suggest  a number  of  nondestructive  testing  applications. 

The  existence  of  thickness  resonances  in  the  reflection  and 
transmission  spectra  of  composites,  analogous  to  those  present  in  monolithic 
materials,  suggest  that  a type  of  generalized  thickness  gauging  can  be 
performed  using  spectrum  analysis.  For  example  the  peak  position  of 
thickness  resonances  below  the  forbidden  band  may  be  given  approximately 
by  the  expression, 

(16)  f «nc/2d  (for  n less  than  the  number  of  laminae) 


where  d is  the  thickness  of  the  composite 
n is  an  integer 

and  c is  the  phase  velocity  of  sound  in  the  composite. 

Clearly,  the  position  of  the  resonance  peak  is  a function  of  the  sample 
thickness  a«d  the  velocity  of  sound,  the  latter  being  a function  of  the 
elastic  constants  (Cjj)  of  the  material  and  its  density. 

For  longitudinal  waves  propagating  normal  to  the  surface  of  a composite 
with  hexagonal  symmetry,  the  velocity  of  sound  in  the  low  frequency  limit 
is  given  by. 


c " 


c33/p 
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Thus  the  resonance  peak  position  In  the  composite  spectrum  can  be  used 
to  monitor  variations  In  density,  thlckh'-ss  and  elastic  constants.  An 
efficient  nwans  of  accomplishing  this  monitoring  it* to  use  the  amplitude 
of  one  component  of  the  frequency  spectrum  at  or  near  a resonance  peak  to 
modulate  the  intensity  of  a c-:>can.  in  this  manner  a composite  sample  may 
be  Imaged  using  a particular  frequency  of  ultrasound.  For  sharp  resonances 
small  changes  In  thickness  or  modulus  will  cause  large  changes  in  the  ampli- 
tude of  ultrasound  being  transmitted  at  or  near  this  frequency,  producing 
corresponding  changes  In  c-scan  Intensity. 


An  example  of  such  a transmission  c-scan  is  shown  in  Figure  9,  which 
is  a c-scan  of  two  boron/aluminum  h ply  unidirectional  tensile  specimens. 

The  nominal  thickness  of  these  specimens  Is  O.Oltl"  and  their  ultrasonic 
thickness  resonances  for  longitudinal  waves  occur  at  about  3.8  HHz. 

(see  Figures  !0a  and  lib).  The  c-scan  Intensity  was  modulated  using  the 
amplitude  of  the  spectral  component  centered  at  3.7  KHz.  In  this  manner 
slightly  thicker  specimens  having  peak  transmission  ai  lower  frequencies 
produce  darker  c-scans  than  thinner  ones.  In  fact,  micrometer  measurements 
of  the  specimens  showed  that  the  darker  specimen  (#1718  Figure  9)  had  an 
average  thickness  of  0.0419",  while  the  lighter  specimen  had  an  average 
thickness  of  0.0402".  Ultrasonic  frequency  spectra  plotted  for  selected 
positions  on  each  specimen  show  that  the  expected  resonance  shifts  do  in 
fact  occur,  the  frequency  shift  between  the  dark  and  light  areas  being 
about  0.2  KHz. 

The  agreement  between  these  two  results  may  be  checked  by  differentiating 
formula  (16)  which  gives  for  the  frequency  shift. 


* %Mk*  ■'  ''“'■’’•W 

This. result  Is  well  within 
ness  measurements. 


0.19  KHz 

the  experimental 


variation  In  the  thick- 


In  addition  to  the  difference  in  average  intensity  observed  between 
the  two  specimens  in  Figure  9 there  Is  a variation  seen  within  each  specimen. 
This  gradation  In  shading  appears  as  bands  oriented  parallel  to  the  fiber 
direction  in  the  specimens.  These  bands  also  have  frequency  shifts 
associated  with  them,  and  the  peak  Intensity  of  the  spectral  thickness 
resonances  associated  with  the  bands  varies  only  slightly.  This  indicates 
that  these  bands  are  associated  with  thickness  or  modulus  variations  rather 
than  delamlnatlons. 

The  sensitivity  of  the  resonances  to  delamlnatlons  can  be  seen  In 
Figure  11,  which  Is  a c-scan  of  a graphIte/epoxy  panel  for  which  the  intensity 
is  modulated  by  the  acoustic  amplitude  corresponding  to  a reflection  thick- 
ness resonance  (ant  I resonance) . This  amplitude  is  small  except  In  the 
presence  of  defects.  Using  this  technique  a number  of  natural  and 
synthetic  defects  were  revealed  without  the  necessity  of  accessing  both 
sides  of  the  sample. 
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A potential  application  of  the  analytical  node!  presented  here  Is 
the  prediction  of  the  frequency  spectra  from  diffusion  bonded  or  adhesively 
bonded  laminates  with  arbitrary  stacking.  Because  the  model  Is  quite 
accurate  for  arrays  of  monolithic  materials,  the  spectrum  resulting  from 
a particular  condition  In  the  laminate  can  be  predicted  without  reference 
to  standard  specimens. 


SUMMARY  AND  CONCLUSIONS 

A model  has  been  constructed  for  describing  the  Interactions  of 
ultrasonic  plane  waves  with  laminated  structures.  Using  this  model 
measurable  parameters  such  as  the  frequency  dependent  transmission  and 
reflection  coefficients  can  be  accurately  calculated  for  laminates 
consisting  of  layers  of  elastic  monolithic  materials. 

Experiments  and  calculations  for  some  simple  finite  periodic 
laminates  have  demonstrated  the  following: 

(1)  In  general  ultrasonic  frequency  spectra  of  laminates  including  periodic 
laminates  are  not  simply  periodic  functions. 

(2)  If  a finite  periodic  laminate  Is  modeled  using  the  model  for  an 
infinite  laminate,  regions  of  high  frequency  attenuation  will  be  localized 
accurately  for  laminates  having  a sufficient  number  of  layers.  However, 
this  model  cannot  predict  detailed  transmission  or  reflection  spectra. 

(3)  For  wavelengths  longer  than  t (where  I is  the  period  of  the 
laminate)  periodic  laminates  behave  acoustically  In  much  the  same  way  as 
dispersive  monolithic  materials  and  their  transmission  and  frequency  spectra 
exhibit  thickness  reasonances,  with  spacing  Af  “ c(f) , where  c(f)  is  the 

2d 

velocity  in  the  composite.  For  wavelengths  on  the  order  of  n1(n  an 
integer)  the  laminate  exhibits  anamalous  dispersion  and  strong  attenuation. 

(4)  Structural  composites  (e.g.  graphIte/epoxy,  flberglas,  B/Al) 
have  acoustic  properties  which  are  qualitatively  similar  to  those  described 
above  for  simple  laminates. 

The  characteristics  of  laminate  spectra  reveal  a number  of  very 
important  facts  which  should  be  considered  when  performing  ultrasonic 
Inspection  of  layered  composite  materials. 

When  inspecting  thick  laminates  ultrasonlcal ly,  the  existence  of 
forbidden  transmission  bands  and  narrow  thickness  resonances  dictates 
that  ultrasonic  frequencies  be  carefully  selected.  This  Is  particularly 
important  If  tuned  transducers  are  being  used.  The  location  of  forbidden 
frequency  bands  will  vary  with  material  properties,  thickness,  lay^-up 
and  the  number  of  plies  in  a laminate  but  will  remain  relatively  constant 
In  thick  laminates  having  the  some  periodicity. 
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Features  of  ultrasonic  spectra  can  be  used  as  sensitive  monitors  of 
changes  in  thickness  and  modulus  both  of  which  are  linearly  related  to 
thickness  resonance  frequencies.  Acoustic  amplitudes  of  waves  propagating 
near  these  resonant  frequencies  can  be  used  to  modulate  c-scan  intensities, 
thereby  generating  maps  of  a sample,  which  reveal  small  changes  in  either 
modulus  or  thickness. 

The  model  which  has  been  constructed  for  uitrasoaic  wave  propagation 
in  laminates  has  proven  to  be  quite  accurate  In  predicting  detailed 
features  of  the  ultrasonic  frequency  spectra  of  layered  structures  and 
should  be  useful  for  describing  wave  propagation  in  any  structure 
consisting  of  monolithic  laminae. 


RECOMMENDATIONS 

1.  The  existing  model  for  wave  propagation  in  laminates  should  be 
applied  to  more  complicated  laminates  including  laminates  with  defects  and 
aperiodic  arrays.  This  modeling  will  require  the  use  of  effective  parameters 
when  dealing  with  viscoelastic  and  fiber  reinforced  materials. 

2.  Efforts  should  be  made  to  generalize  the  present  model  to  cases 
involving  scattering,  shear  wave  propagation  and  propagation  of  waves  at 
oblique  incidence  to  laminate  interfaces. 

3.  Ultrasonic  NOT  techniques  and  devices  should  be  developed  which 
exploit  the  resonant  nature  of  laminated  structures.  These  might  include 
the  following: 

a.  Selection  of  frequencies  for  maximum  penetration  Into  laminates. 

b.  The  use  of  tuned  transducers,  filters  and  amplifiers  for  Improved 
signai-to'noise  ratio  when  interrogating  laminated  composites  near  their 

resonant  frequencies. 

c.  The  use  of  mult  I -frequency  techniques  for  monitoring  lay-ups 
and  consistency  in  structural  composites. 
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VIRTUAL  INTERFACE 
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Figure  2 Ultrasonic  Pulse  Incident 


upon  Layered  Media 


Frequency  In  HH* 

Transmission  Spectrum  for  0.1 
Steel  Plate  In  Water 


Frequency  In  MH* 

Refloctfon  Spectrum  for  O.^n"  Thick 
Steel  Plate  In  Water 


Transmission  Spectrum  for  two  1.22mm  Glass  Plates  In  Water 
(0.1mm  separation  ) 


Plates  In  Water 
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Figure  k Apparatus  for  Measuring  Acoustic  Reflection 
and  Transmission  Coefficients  as  a Function 
of  Frequency 
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(Measured  Experimentally) 


Ultrasonic  Frequency  Spectra  of  Graphite  Fpoxy  Panels 


Fftjure  8 Ultrasonic  Velocity  vs  Frequency  I 


rasonJc  Frequency  Snectra  for  Roron  Aluminum 
(letters  refer  to  Flqure  9) 


Figure  lOb  Ultrasonic  Frequency  Spectra  for  Boron  Aluminum 
(letters  refer  to  Figure  9) 


Figure  11  C-scan  for  Graphite  Fooxy  Panel  Showino 
Natural  and  Synthetic  Flaws 
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